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Abstract 

Following the previous letter on the first microscopic description of the antimagnetic rotation 

O 

(AMR) in 105 Cd, a systematic investigation and detailed analysis for the AMR band in the frame- 



work of tilted axis cranking (TAC) model based on covariant density functional theory are carried 



J> | out. After performing the microscopic and self-consistent TAC calculations with an given density 

' functional, the configuration for the observed AMR band in 105 Cd is obtained from the single- 

00 ' 

particle Routhians. With the configuration thus obtained, the tilt angle (fo for a given rotational 

£^ , frequency is determined self-consistently by minimizing the total Routhian with respect to the an- 

(N ■ 

gle 0q. In such a way, the energy spectrum, total angular momenta, kinetic and dynamic moments 
of inertia, and the B(E2) values for the AMR band in 105 Cd are calculated. Good agreement with 

X - 

the data is found. By investigating microscopically the contributions from neutrons and protons to 
the total angular momentum, the "two-shears-like" mechanism in the AMR band is clearly illus- 
trated. Finally, the currents leading to time-odd mean fields in the Dirac equation are presented 
and discussed in detail. It is found that they are essentially determined by the valence particles 
and/or holes. Their spatial distribution and size depend on the specific single-particle orbitals and 
the rotational frequency. 

PACS numbers: 21.60.Jz, 21.10.Re, 23.20.-g, 27.60. +j 
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I. INTRODUCTION 

In the previous letter jlj], we reported the first fully self-consistent and microscopic investi- 
gation of antimagnetic rotation (AMR) in 105 Cd in the framework of the tilted axis cranking 
(TAC) model based on the covariant density functional theory (CDFT). We found that the 
experimental data of this band 2J have been reproduced rather well without any adjustable 
parameters. The present paper goes in more details along this direction by providing a series 
of interesting results for the AMR band in 105 Cd and the microscopic interpretation of AMR 
in nuclei. 

In general, rotational bands are built on intrinsic states violating rotational symmetry. 
In the common cases this symmetry is broken by a substantial quadrupole deformation and 
therefore these bands show strong electric quadrupole {E2) transitions between the rota- 
tional states. Many exciting phenomena have been discovered and predicted in this field such 
as backbending 3], alignment phenomena j^, 5], and superdeformed rotational bands jsj. 
Such bands are most commonly interpreted as arising from the coherent collective rotation of 
many nucleons around an axis perpendicular to the symmetry axis of the deformed density 
distribution Q. 

Over the past decades, however, a new type of rotational band with strong magnetic dipole 
(Ml) and very weak E2 transitions has been discovered experimentally in nearly spherical 
light lead isotopes and other groups of nuclei (for reviews see Refs. 8l-ll0l|). The intriguing 
feature here is that the orientation of the rotors is not specified by the deformation of the 
overall density but rather by the current distribution induced by specific nucleons moving 
in high-j orbitals. 

The explanation of such bands in terms of the "shears mechanism" was firstly given 
in Ref. where a decreasing tendency of B(M1) values with total angular momentum 
was also predicted. In this interpretation, the angular momentum vectors of the high-j 
proton and neutron form two blades of a pair of shears and are almost perpendicular to each 
other at the bandhead. Along the bands, energy and angular momentum are increased by 
closing the blades of the shears, i.e., by aligning the proton and neutron angular momenta. 
Consequently, rotational bands are formed in spite of the fact that the shapes of these nuclei 
stay nearly spherical. A clear evidence for this new rotation mode has firstly been provided 
through lifetime measurements for four Ml bands in 198 > 199 Pb 12J. In order to distinguish 
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this kind of rotation from the usual collective rotation in well-deformed nuclei (called electric 
rotation), the name "magnetic rotation" 13] was introduced. It alludes to the fact that here 



the magnetic moment is the order parameter inducing a violation of rotational symmetry 
and thus causing rotational-like structures in the spectrum [14 1. This forms an analogy to 
a ferromagnet where the total magnetic moment, the sum of the atomic dipole moments, is 
the order parameter. 

In an antiferromagnet, on the other hand, one-half of the atomic dipole moments are 
aligned on one sublattice and the other half are aligned in the opposite direction on the 
other sublattice. In such a way, the net magnetic moment in an antiferromagnet is canceled 
out. However, it is still an ordered state since the isotropy of such a state is also broken like 
a ferromagnet. 

In analogy with an antiferromagnet, a similar phenomenon can be predicted in nuclei for 
"antimagnetic rotation" (AMR) in Ref. jsilllj]: in specific nearly spherical nuclei, subsystems 
of valence protons (neutrons) are aligned back to back in opposite directions and nearly 
perpendicular to the orientation of the total spin of the valence neutrons (protons) (see 
Fig. IE]). Such an arrangement of the proton and neutron angular momenta also breaks 
rotational symmetry in these nearly spherical nuclei and causes excitations with rotational 
character on top of this bandhead. Along this band, energy and angular momentum are 
increased by simultaneous closing of the two blades of protons and neutrons toward the 
total angular momentum vector. Consequently, a new kind of rotational bands in nearly 
spherical nuclei is found showing some analogy with an antiferromagnet in solid state physics. 

AMR is expected to be observed in the same regions in the nuclear chart as magnetic 
rotation 9j. However, it differs from magnetic rotation in two aspects. First, there are no Ml 
transitions in the AMR band since the transverse magnetic moments of the two subsystems 
are antialigned and canceled out against each other. Second, as the antimagnetic rotor is 
symmetric with respect to a rotation by 180° about the angular momentum axis, the energy 
levels in the AMR band differ in spin by 2h and are connected by weak E2 transitions 
reflecting the nearly spherical core. Moreover, the phenomenon of AMR is characterized 
by a decrease of the B(E2) values with spin, which has been demonstrated by lifetime 
measurements 16J • Since AMR was proposed [9| , it has attracted more and more attention 



both experimentally and theoretically. To^ate, experimental evidences of AMR have been 



reported in Cd isotopes including 105 Cd [2|, 106 Cd [16|, 108 Cd (13,1131, and 110 Cd [19|. In 
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addition, the occurrence of this phenomenon in 109 Cd 20], 100 Pd 21], 144 Dy 22] still needs 
further investigations. 



Theoretically, the AMR bands have been discussed in simple geometry 



as the TAC model 



0,Q : 



17 



20 



as well 



2l| . The TAC model can explicitly construct classical 



vector diagrams showing the angular momentum composition. This is of great help in 
visualizing the structure of the rotational bands. The quality of the cranking approximation 



for principal axis cranking 



23( | . tilted axis cranking 24J, and aplanar tilted axis cranking 25] 



has been discussed and tested within the particle rotor model. Based on the TAC model, 
many applications have been carried out in the framework of the pairing plus quadrupole 



model 



20] or the microscopic-macroscopic model 



16 



17 



2l|. 



In these investigations, however, the polarization effects which are expected to have a 
strong influence on the quadrupole moments and thus the B(E2) values are either neglected 
completely or taken into account only partially by minimizing the rotating energy surface 
with respect to a few deformation parameters. Moreover, the nuclear currents, which are 
the origin of symmetry violation in nuclei with AMR, are not treated in a self-consistent 
way in these models. Therefore, it is evident that a full understanding of AMR requires 
self-consistent and microscopic investigations including all degrees of freedom and based on 
reliable theories without additional parameters. Such calculations are not simple, but they 
are nowadays feasible in the framework of density functional theories. 

During the past decades, the CDFT has received wide attention due to its success in 
describing many nuclear phenomena in stable as well as in exotic nuclei 2614281] . On the basis 
of universal functionals and without any additional parameters, the rotational excitations 
can be described in practical applications within the self-consistent cranking relativistic 
mean-field (RMF) framework. In rotating nuclei, time reversal symmetry is broken by the 
Coriolis operator and by unpaired nucleons. This leads to time-odd components in the fields. 
In non-relativistic density functionals, these time-odd components are usually not very well 
defined, because the phenomenological parameters of these functionals are determined by 
experimental ground state properties in even-even nuclei where these time-odd fields vanish. 
In covariant density functionals this problem does not occur. The time-odd fields in the 
Dirac equation correspond to non-vanishing space-like components of the vector potential. 
They are often called nuclear magnetism and are induced by the various currents. Lorentz 
invariance requires that the time-like and the space-like components of the vector fields are 
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determined by the same coupling constants, which are well determined by experimental 



ground state properties. Dobaczewski et al 



29[ have emphasized that Galileian invariance 



induces constraints for the coupling constants in non-relativistic theories too. However, 
in many of the successful parameterizations of non-relativistic density functionals these 
conditions are not fulfilled and the spin parts of the functionals which are also connected 
with time-reversal breaking are not influenced by Galileian invariance. Therefore, it seems 
to us very important to use relativistic functionals for the description of magnetic and 
antimagnetic rotational bands, where the symmetry breaking mechanism is introduced by 
the currents. 

A code for the cranked RMF model with arbitrary orientation of the rotational axis, 
i.e., three-dimensional cranking, has been developed in Ref. 30] . However because of the 
numerical complexity of this code, so far, this model has been applied only for the magnetic 
rotation in 84 Rb. Focusing on more investigations of bands with magnetic rotation, a com- 
pletely new computer code for self-consistent two-dimensional cranked RMF theory has been 
developed in Ref. 31]. It includes significant improvements in implanting the simplex sym- 
metry V y T, quantum number transformation, and orientation constraints. Very recently, 
the TAC model based on a relativistic point-coupling Lagrangian which allows considerable 
simplification was established and applied successfully to the mag netic rotation both in light 

60 T 



nuclei such as eu Ni 32] and in heavy nuclei such as 198 > 199 Pb 33]. 

The first fully self-consistent and microscopic investigation of AMR has been reported 
recently in the framework of TAC based on a relativistic point-coupling Lagrangian in CDFT 
in Ref. The present paper focuses on more detailed results in 105 Cd including the 
Routhians, energies, angular momenta, moments of inertia, electromagnetic transitions, as 
well as the nuclear currents and densities. Moreover, the microscopic interpretation of 
AMR in nuclei such as the configuration, the characteristics of AMR, the "two-shears-like" 
mechanism, and the alignment of the angular momenta are discussed in detail. In Sec. 
II, the TAC formalism is implemented in CDFT using a point-coupling Lagrangian. The 
numerical details are given in Sec. III. In Sec. IV, a series of calculated results of the AMR 
band in 105 Cd are presented. Finally, a summary is given in Sec. V. 
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II. THEORETICAL FRAMEWORK 



The basic building blocks of point-coupling vertices are 

(tpOTip), Oe{l,r}, re{l, TA4 ,7B,757M><v}> (!) 

where ip is the Dirac spinor field of nucleon, f is the vector of isospin Pauli matrices, and 
T generally denotes 4x4 Dirac matrices. Throughout this paper, vectors in isospin space 
are denoted by arrows and space vectors by bold type. Greek indices /i and v run over the 
Minkowski indices 0, 1, 2, and 3. 

A general effective Lagrangian can be written as a power series in ipOTip and their 
derivatives, with higher-order terms representing in-medium many-body correlations. In 
the actual application we start with the following Lagrangian density of the form: 

C = C bce + £ 4f + £ hot + £ dcr + £ em , (2) 
= ip^^d^ — m)ip 

-iF^F^-el-^WMv (3) 

which includes the Lagrangian density for free nucleons £ free , the four-fermion point-coupling 
terms £ 4f , the higher order terms £ hot accounting for the medium effects, the derivative terms 
£ der to simulate the effects of finite-range which are crucial for a quantitative description for 
nuclear density distributions (e.g., nuclear radii), and the electromagnetic interaction terms 
£ em . The higher order terms lead in the mean field approximation to density dependent 
coupling constants with a de nsity dependence of polynomial form. Lagrangians with more 



general density dependencies 34j can be used in a similar way. 

In the tilted axis cranking approach based on CDFT, it is assumed that the nucleus 
rotates around an axis in the xz plane and the Lagrangian in Eq. (j2J) is transformed into a 
frame rotating uniformly with a constant rotational frequency, 

0=^,0,^) = (f2cos0 n ,O,nsin0 n ), (4) 



6 



where 6q := <(fl, e x ) is the tilt angle between the cranking axis and the x axis. From this 
rotating Lagrangian, the equation of motion for the nucleons can be derived in the same 



manner as in the meson-exchange case 35|, |36| and one finds 



[a • (-zV - V) + P(m + S) + V° - SI ■ J]^ fc = e k i/> k , (5) 

where J = L + |S is the total angular momentum of the nucleon spinors, and represents 
the single-particle Routhians for nucleons. The relativistic fields S(r) and V^(r) read 

S(r) = a s p s + PsPs + IsPs + &s&Ps, (6a) 
V°(r) = a v pv + IvPv + &v^Pv + T 3 a T vPTV + t 3 5 T v^Ptv + eA°, (6b) 
V(r) = a y jy + iv(jv) 3 + Sv^jv + T 3 a TV j TV + r 3 5 TV Aj TV + eA, (6c) 

with e the electric charge unit vanishing for neutrons. As usual, it is assumed that the 
nucleon single-particle states do not mix isospin, i.e., the single-particle states are eigenstates 
of r 3 . Therefore only the third component of isovector potentials survives. The Coulomb 
field A°(r) is determined by Poisson's equation 

- AA°(r) = e Pc . (7) 



From prior experience in the meson-exchange case [35l . |37J , the spatial components of the 
electromagnetic vector potential A(r) are neglected since these contributions are extremely 
small. 

In comparison with the coupled equations of motion in the case of meson-exchange po- 
tentials, one finds in the point-coupling model considerable simplifications. In particular, 
apart from the Laplace equation ((7j), one does not have to solve the Klein-Gordon equations 
for mesons. 

Since the Coriolis term f2 • J in the Dirac equation (JHJ) breaks time reversal symmetry 
in the intrinsic frame, currents are induced and as a consequence spatial components of the 
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vector potential V(r). The densities and currents in Eqs. ([6]) have the form 

A 

^^(r)^(r), 

i=l 

^^i(r)^i(r), 

i=i 

^^t( r )a^(r), 



py(rj 



i=l 
4 



i=l 
i=l 

X>1 



r)T 3 ^»(r), 

r)ar 3 ^(r 
^l-r 3 



(8a) 
(8b) 
(8c) 
(8d) 
(8e) 
(8f) 



i=l 



Following, as usual, no-sea approximation the sums are taken over the particles states in 
the Fermi sea only, i.e., the contribution of the negative-energy states are neglected. In fact 
these contributions resulting in vacuum polarization are not neglected completely because 
this effect is taken to account in a global way by the adjustment of the parameters in the 
Lagrangian to experimental data 



By solving the equation of motion iteratively, one finally obtains the total energy in the 
laboratory frame 



Etot = E kin + E- mt + E t 
which is composed of a kinetic part 



\~ E c .m. i 



(9) 



an interaction part 

Eint 



A 

E kin = d^rY^ W [a-p + Pmtyi 

i=l 



+ ^ a v(pv - 3 ■ j) + \^Tv{p 2 TV ~ hv ■ 3tv) 
+\lv(p 2 v - 3 ■ 3? + ^S v (p v Ap v - jAj) 



(10) 



+ ^ TV (p TV Ap TV - jrvAjrv) f , 



(11) 



an electromagnetic part 



' d 3 r^eA pc, (12) 



and the center-of-mass (cm.) correction energy £ c . m , accounting for the treatment of center- 
of-mass motion 

= -^X( P L), (13) 

where A is the mass number and P c . m . = Ylt Pi 1S t ne total momentum in the center-of-mass 
frame. It has been shown that this cm. correction based on the self-consistently determined 
wave functions provides more reasonable and reliable results than that based on the simple 



oscillator model 



For each rotational frequency Q the expectation values of the angular momentum com- 
ponents J = (J x , J y , J z ) in the intrinsic frame are given by 

A 

8=1 

Jy = 0, (14b) 

A 

J* = <£> = £rf\ (14c) 
i=i 

and by means of the semiclassical cranking condition 

J = \/(l) 2 + (I) 2 = vW^l), (15) 

one can the relate absolute value of the rotational frequency Q to the angular momentum 
quantum number / in the rotational band. 

The orientation of the angular momentum vector J is represented by the angle 8j : = 
<( J, e x ) between the angular momentum vector J and the x axis. In a fully self-consistent 
calculation, the orientation 6j of the angular momentum J should be identical to the ori- 
entation 6q of the angular velocity ft. 

The quadrupole moments Q 2 q and Q 2 2 are calculated by 



(16a) 



^(^-y 2 ), (16b) 
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and the deformation parameters (3 and 7 can thus be extracted from 

(3 = /c 



1 



arctan 



In 1 

M 22i 
Q20 



by using the relations 



3A 2 

Q20 — ~^RqO,20) 



Q 



22 



3A 
Air 



RqO.22, 



(17a) 
(17b) 

(18a) 
(18b) 



with R = 1.2 A 1 / 3 fm. Note that the sign convention in Ref. [42| is adopted for the definition 
of 7 here. 



The nuclear magnetic moment is, in units of the nuclear magneton, given by 43] 



ft 



A 

Wd»r 
i=i J 



mc 



he 



#J(r)r x ati[)i(r) + Kip] (r) /3E^i(r) 



(19) 



where the charge q (q p = 1 for protons and q n = for neutrons) is given in units of e, m 
the nucleon mass, and k the free anomalous gyromagnetic ration of the nucleon (n p = 1.793 
and n n = —1.913). 

From the quadrupole moments and the magnetic moment, the B(M1) and B(E2) tran- 
sition probabilities can be derived in semiclassical approximation 

\2 



3 3 

B{M1) = —nl = —(fi x sm6j-fi z cos6jY 

07T 07T 



B(E2) 



Q P 2Q cos 2 9 j 



Q p 22 (l + sm 2 6j) 



(20a) 
(20b) 



where Q?, and Q22 corresponds to the quadrupole moments of protons. 

As in the tilted axis cranking relativistic mean-field (TAC-RMF) model with meson- 
exchange interaction 3l|, the directions of the principal axes of the density distribution are 
fixed along the x, y, and z axes by a quadratic constraint on (Q2-1) = in order to eliminate 
the numerical instabilities which results from the freedom of rotations around the y axis. 
More details on this treatment can be seen in Ref. 31 1. 



III. NUMERICAL DETAILS 



In the present work, as in Ref. 3l|, the invariance of space reflection V and the combi- 
nation of the time reversal and reflection in the y direction V y T are used in the code for 
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the solution of the TAC-RMF model with a point-coupling interaction. The Dirac equation 
for the nucleons is solved in an isotropic three-dimensional Cartesian harmonic oscillator 
basis with Nf = 10 major shells and the oscillator frequency given by hj = 41A" 1 / 3 . By 
increasing Nf from 10 to 12, the changes of total energies and total quadrupole moments 
are within 0.1% and 4% for the ground state of the nucleus 105 Cd, respectively. Therefore, 
a basis of 10 major oscillator shells is adopted in the present calculations. Moreover, the 



point-coupling interaction PC-PK1 [44j is used for the Lagrangian without any additional 
parameters and pairing correlations are neglected. 

In the recent experiment of 105 Cd a negative-parity band was reported in forming the 
yrast line above I = 23 /2h. It was assumed that this band corresponds to three aligned 
neutrons with the configuration ^[^iiMfi^) 2 ] cou Pl e d to a pair of g Q / 2 proton holes. As 
in Ref. j^j, in the present calculations the configuration is fixed at the bandhead: The odd 
neutron is kept in the lowest level of the hu/2 shell. The remaining nucleons are treated self- 
consistently by filling the orbitals according to their energy from the bottom of potential. 
This automatically leads to the configuration shown in Fig. [TJ 

In order to keep the configurations unchanged with increasing rotational frequency Q, 



we adopt the same prescription as in Ref. 



311 ] : starting from the blocked level \ipi(Q)) 



at frequency f2, we block for frequency Q + 8Q the orbit j which maximizes the overlap 

fain + styMn)), i.e., 

(tl) j (n + 5n)\t(>i(n)) = i + o(5n). (21) 

In the practical calculations, the rotational frequency Q is increased successively in steps of 
5(hQ) = 0.1 MeV and in each case the prescription Eq. (J2T1) is applied for the occupation of 
the various levels. Of course, as we see from Fig. [TJ this is in our case equivalent to apply, 
for each Q, this prescription only for the odd neutron in the hu/2 shell and to determine the 
occupation of the other levels in a self-consistent way by filling the levels from the bottom 
of the well. 



IV. RESULTS AND DISCUSSION 

In this section we present the results of the TAC-RMF calculations for the nucleus 105 Cd, 
in particular the Routhians, energies, angular momenta, moments of inertia, electromagnetic 
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transitions, the currents and densities. Moreover, the microscopic interpretation of AMR in 
nuclei such as the configuration, the characteristics of AMR, the two-shears-like mechanism, 
and the alignment of the angular momenta are discussed in detail. 



A. Routhians 
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FIG. 1: (Color online) Single neutron (upper panel) and single proton (lower panel) Routhians 
near the Fermi surface in 105 Cd as functions of the rotational frequency. The filled circles indicate 
the occupied levels, while the open circles in the lower panel represent the holes in the g$/2 shell. 

At finite cranking frequency Q, the Coriolis term violates time reversal invariance in the 
intrinsic frame and this leads to an energy splitting of the degenerate time reversal conjugate 
states at Q = 0. This can be seen in Fig. [1] where the single particle Routhians e« in Eq. (JHJ) 
near the Fermi surface are shown as functions of the rotational frequency. The amplitudes 
of this energy splitting in 105 Cd range from to 2.3 MeV. As expected, we find a larger 
splitting for states with larger expectation values \j x \. This results from the Coriolis term 
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which reads Q X J X when the rotational axis is along the x axis. It should be noted that even 
without Coriolis term, i.e. for zero angular velocity, time reversal symmetry is broken by the 
current of the odd particle. Therefore the energy levels are not exactly twofold degenerate 
and some of them show some splitting in the energy already at Q = 0. This phenomenon 
has been discussed, e.g., in Ref. [43]. 

In Fig. [TJ the occupied single-particle levels are marked by filled circles. For the protons, 
one can distinctly see that two proton holes are always sitting at the top of the gg/2 shell. For 
the neutrons, there are totally seven valence neutrons above the N = 50 shell. As discussed 
in the previous section, the odd neutron is kept fixed in the lowest level in the fan/2 shell. The 
other six remaining neutrons are treated self-consistently by filling the orbitals according to 
their energy and distributed over the (g? 72^5/2) shell. We find that there is a strong mixing 
between these orbitals. In this sense, we have here, in the present microscopic calculations, 
a similar configuration as that of Ref. 
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FIG. 2: (Color online) TAC-RMF results for total Routhians (solid lines) as functions of the tilt 
angle 0q. The rotational frequencies range from Ml = 0.2 MeV to 0.7 MeV. The solid dots represent 
the minima for the corresponding frequency. 

So far we kept the tilt angle 6q always fixed at 9q = 0. Of course, this angle should be 
determined in a self-consistent way by minimizing the total Routhian 



tot 



(cos 6 n QJ x + sin 6 n QJ z ) 



(22) 



with respect to the angle 9q. Therefore we repeated this calculations for a range of 8^ values 
and in Fig. |2] we show the total Routhians as functions of the tilt angle #n for rotational 
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frequencies ranging from hVt = 0.2 MeV to 0.7 MeV. It turns out that the minimum of the 
total Routhians (solid dot) remains for all angular velocities Q at the tilt angle 8q = 0°. 
This demonstrates clearly that the rotational axis is always parallel to the x axis in the 
present antimagnetic rotor. Of course this is a natural consequence of the symmetry for 
the antimagnetic configuration, but it is not trivial, because symmetries can be broken in 
the self-consistent solution. Of course, having this result we can conclude that, in principle, 
both the principal axis cranking (PAC) and TAC calculations can be used to describe AMR. 
The reasons for carrying out TAC calculations here is that the relative orientation of the 
g 9 / 2 proton hole vectors can be calculated easily (see Fig. [7] below). 



B. Energy and moment of inertia 
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FIG. 3: (Color online) Energy spectrum obtained from the TAC-RMF calculations in comparison 
with the data (upper panel) as well as the corresponding contributions from kinetic energy, inter- 
action energy, Coulomb energy, and center-of-mass correction energy (lower panel). These results 
are normalized at / = 23 /2h. 



In the upper panel of Fig. [3l the calculated total energy (full line) is shown in comparison 
with the data 2j. The experimental energy spectrum is reproduced in an excellent way 
by the present self-consistent calculations. In Eq. ((9]), the calculated energy is divided into 
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kinetic energy, interaction energy, Coulomb energy, and center-of-mass correction energy. 
The contributions to the total energy from these four parts are shown in the lower panel of 
Fig. El It is found that the Coulomb energy and center-of-mass correction energy are barely 
changed in amplitude (less than 0.5 MeV). In comparison, the kinetic energy increases 
monotonously up to roughly 6 MeV and the interaction energy declines monotonously up to 
roughly 2 MeV. Therefore we can conclude that the rotational excitation energy in the band 
comes mainly from the kinetic energy with a moderate competition from the interaction 
energy. 
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FIG. 4: (Color online) The total angular momentum (upper panel) as well as the kinetic 3^ 
and dynamic 0^ moments of inertia (lower panel) as functions of the rotational frequency in the 
TAC-RMF calculations in comparison with the data. 



In Fig. HI the total angular momentum (upper panel) as well as the kinetic 3 and 
dynamic 3^ moments of inertia (lower panel) are shown as functions of the rotational 
frequency and compared with the data. Taking into account the quantal corrections 24], 
the total angular momentum J calculated in TAC model corresponds to the quantum number 
of the angular momentum 1+1/2 since + 1) ~ 1 + 1/2. This prescription permits us to 

compare the TAC calculations with the experimental data. In the present AMR case, states 
differing by two units of angular momentum are arranged into a AI = 2 band. Therefore, 
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the experimental rotational frequency is extracted as 

J = I- 1/2, nfiexp = ^ 7 (/-H-2). (23) 

Meanwhile, the experimental values of kinetic ^ and dynamic moments of inertia 
are calculated from the transition energies by using the finite difference approximation for 
AI = 2 bands, 

^ = E 1 (I + 2^I)-E 1 (I ^1 -2Y (24b) 
with nn cxp = \E^{I ->■ I — 2). 

In the upper panel of Fig. [4j it can be found that the calculated total angular momenta 
increase almost linearly with increasing rotational frequency and agree with the data very 
well. This leads to fairly constant values for both the kinetic and dynamic moments of inertia 
as shown in the lower panel. Generally, the fact that both the kinetic and the dynamic 
moments of inertia stay roughly constant is well reproduced. However, the calculated values 
for 3^ underestimate the data. 



C. Transition probability B(E2) and deformation 

Typical characteristics of AMR include weak E2 transitions, reflecting the small defor- 
mation of the core, which results in large ratios of the dynamic moments of inertia 0^ to 
the reduced transition probability B(E2) values. What's more, the B(E2) values rapidly 
decrease with the angular momentum. In the upper and lower panels of Fig. |5l the calculated 
B(E2) values and the 3 (2) / 1 B(E2) ratios are shown as functions of the rotational frequency 
in comparison with the data. It is found that the B(E2) values are in very good and the 
Z^/B(E2) ratios in reasonable agreement with the data. Quantitatively, one can see that 
the calculated B(E2) values are small (< 0.14 e 2 b 2 ) and this leads to the large 3^ /B(E2) 
ratios (> 150 MeV _1 ^. 2 e~ 2 &~ 2 ). Moreover, the fact that the B(E2) values decrease smoothly 
with the growing rotational frequency, which results in an increase of the Z {2) /B(E2) ratios, 
is connected with the interpretation of the two-shears-like mechanism. 

It should be noted that in some nuclei (e.g., 108 Sn and 109 Sb) in the A ~ 110 mass 
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451 ] . a so-called "smoothly terminating band" is observed. Its structure gradually 
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FIG. 5: (Color online) B(E2) values (upper panel) and 3^ / B(E2) ratios (lower panel) as functions 
of the rotational frequency in the TAC-RMF calculation in comparison with the data. 

develops from a prolate shape rotating around an axis perpendicular to the symmetry axis, 
through a sequence of triaxial shapes, to an oblate shape with the symmetry axis parallel to 
the total angular momentum vector. In this sense, the band terminates when the angular 
momenta of the valence particles are fully aligned, just as the band of the antimagnetic 
rotor terminates when the shears blades are aligned. However, one of the most important 
feature for AMR is the rise of Z^/B(E2) ratios reflecting the fact that as shown in 
the lower panel of Fig. HJ is essentially constant whereas the B(E2) values decline along the 
band. This is not the case for a smoothly terminating band which has an almost constant 
/B(E2) ratio, indicating that the 3^ have falling behavior as similar as B(E2) with the 
increasing spin. 

The decrease of the B(E2) values can be understood by the changes in the nuclear 
deformation. In Fig.[6]we present the evolution of the deformation for neutrons, protons and 
the whole nucleus in the (/?, 7) plane with increasing rotational frequency. It is found that 
the deformations for neutrons, protons and the whole nucleus behave in a similar way, i.e., a 
rapid decrease in j3 together with a small and nearly constant triaxiality. In particular, the 
calculated proton deformation reflecting the charge distribution and therefore responsible 
for the B{E2) values decreases from /3 = 0.18 to /3 — 0.12 with rather small triaxiality 
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p 

FIG. 6: (Color online) Evolutions of the neutron, proton and total deformations in the 7) 
plane with rotational frequency in the TAC-RMF calculations. The arrow indicates the increasing 
direction of rotational frequency. 

(7 < 7°). This is related to the decrease of the B(E2) values shown in the upper panel of 
Fig. [5] and one can thus conclude that the alignment of the shears blades in AMR, i.e., the 
two-shears-like mechanism, is accompanied by a transition from a prolate towards a nearly 
spherical shape. 

D. Two-shears-like mechanism 

In order to examine the two-shears-like mechanism for the AMR band in 105 Cd, we show 
in Fig. [7] the angular momentum vectors of the neutrons J v and the two gg/2 proton- holes 
j n for rotational frequencies ranging from 0.3 to 0.7 MeV. The neutron angular momentum 
J v is defined as 

N N 

J„HX) = ^ n) =5>|J|n>, (25) 

71=1 n=l 

where the sum runs over all the neutron levels occupied in the cranking wave function in 
the intrinsic system. At the bandhead (Ml = 0.3 MeV), each of the two proton angular 
momentum vectors j n are pointing opposite to each other and are nearly perpendicular 
to the vector J v . They form the blades of the two shears. As the rotational frequency 
increases, the gradual alignment of the vectors j n of the two gg/2 proton- holes toward the 
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FIG. 7: (Color online) Angular momentum vectors of neutrons J u and the two 59/2 proton holes 
j w at rotational frequencies from 0.3 to 0.7 MeV. 

vector J v generates angular momentum, while the direction of the total angular momentum 
stays unchanged. This leads to the closing of the two shears simultaneously by moving one 
blade towards the other. In such a way, the two-shears-like mechanism in 105 Cd is clearly 
demonstrated in Fig. [7J Similar results on the distribution of angular momenta for the 



#9/2 proton holes and the /in/2 neutron(s) in other Cd isotopes can a 



so be obtained with 



20| or the microscopic- 



TAC in the framework of either the pairing plus quadrupole model 
macroscopic model Q • 

In a fully self-consistent and microscopic calculation, both the energy and angular momen- 
tum come from individual nucleons. In order to investigate this two-shear-like mechanism in 
a more microscopic way, we show in Fig. [S] the angular momentum alignment for the seven 
neutron particles [panel (a)] and two proton holes [panel (c)] as well as the compositions of 
the neutron [panel (b)] and proton [panel (d)] angular momenta. 

It is found that the spin contributions from the proton and neutron cores, i.e., from the 
orbitals below Z = 50 and A" = 50, are quite small. This 100 Sn core stays relatively inert. 
Instead, the contributions to the angular momenta mainly arise from the high-j orbitals, 
i.e., from proton holes in the c/9/2 shell as well as from neutron particles in the hu/ 2 and g 7 / 2 
shells. 
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FIG. 8: (Color online) Angular momentum alignment for the seven neutron particles [panel (a)] 
and two proton holes [panel (c)] as well as the compositions of the neutron [panel (b)] and proton 
[panel (d)] angular momenta. 

For the protons, the angular momentum is mainly built from the two holes at the top of 
the c/9/2 shell while the spin contribution from the Z = 50 core is rather small (< 0.5h). As 
shown in Fig. U\ the two proton holes cancel each other in the z direction. With growing 
angular velocity, the two proton holes align gradually along the rotational axis and the 
angular momentum in x direction increases. The blades of the two shears are closing. 

For the neutrons, the angular momentum along the x axis originates mainly from the 
neutron particles filling in the bottom of the /in/2 and #7/2 shell whereas the contribution 
from the N = 50 core is less than lh. A neutron sitting in the hu/2 orbit contributes an 
angular momentum of roughly 5H along the x axis. As the rotational frequency increases, the 
contribution of this neutron does not change much and the increase of angular momentum 
is mostly generated by the alignment of the other six neutrons in the g 7 / 2 and d 5 / 2 orbitals. 
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The small bump shown in the spin contributions of the neutron core at Ml ~ 0.5 MeV is 
associated with a virtual crossing between the core g$/2 and the valence (#7/2^5/2) orbitals 
where the spin contributions can not be split properly between these crossing orbitals. As 
the level crossing involves only the filled orbitals, there is no influence on the total angu- 
lar momentum but we see it only in the somewhat artificial division of spin contributions 
between the core and valence particles. 
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FIG. 9: (Color online) Probability variations as functions of the rotational frequency for the main 
components of the valence single-particle states at Ml = 0.7 MeV. The probabilities of j x for seven 
valence neutron orbitals and the probabilities of j z for two valence proton orbitals are displayed in 
the left and right panels respectively. 



In the present microscopic calculation, when a configuration is labeled as having n par- 
ticles, e.g., in the (57/2^5/2) subshells, it does not mean that the particles in the deformed 
rotating potential are in these pure sub-shells, but instead they are distributed over spherical 
orbits with the quantum numbers {n£jj x ), where n is the radial quantum number, (£j) are 
the angular momentum quantum numbers and j x is the magnetic quantum number along 
the x-axis. Therefore, both the alignment of proton holes and neutron particles arise from 
a considerable mixing between the single-particle orbitals. 

In order to quantify these statements somewhat we expand the wavefunction of the ro- 
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tating particle i in terms of the spherical oscillator basis states 

\i)=^2(nijj x \^\n£jj x ) (26) 

and show in Fig. [9]the probabilities \(n£jj x \i)\ 2 that the orbit (n£jj x ) (labeled by the letters 
A, B, C, D, E, F, G) is occupied by a particle i in the following way: In the left panel we 
concentrate on the seven neutrons orbits in the valence shell. Starting with the maximally 
aligned situation at the angular velocity Ml = 0.7 MeV we calculate for each orbit {n£jj x ) 
the particle % with the highest probability to sit in this shell. It turns out that the aligned 
l/in/2, jx = 11/2 orbit is occupied by the odd neutron with a very high probability of roughly 
85 %. In the rest of the valence shell there is more mixing. Of course, these probabilities 
depend on the angular velocity. Next we go to lower frequencies and determine for each 
orbit (n£jj x ) the particle i with the maximal overlap to the corresponding particle at the 
neighboring next higher frequency applying the rule given in Eq. (12T]) and we find that there 
is not much change in the occupation of the aligned lhn/2,j x = 11/2 orbit. At the bandhead 
it is still occupied with a probability of 77 %. As shown in Fig. [8] it always contributes ~ 5H 
to the angular momentum along the x axis. However, the other probabilities decrease with 
decreasing angular velocity. When we reach the bandhead, the corresponding probabilities 
have decreased by more than 20 percentage points. Starting from the bandhead, the increase 
of angular momentum is generated mostly by the alignment of the six neutrons in the gj/2 
and c?5/2 orbit als. 

Next we consider the protons in the right panel of Fig. [91 Here we show only the two 
lgg/2 orbits with the angular momentum quantum numbers j z = ±9/2 along the z-axis. 
both probabilities are equal because of the symmetry with respect to a rotation of 180° 
around the x-axis. At the band head, at Ml = 0.3 MeV, the two protons holes occupy these 
orbitals with 86 % probability, but with increasing angular velocity this occupation rapidly 
decreases because the holes align along the rotational x-axis. This alignment is characterized 
by a considerable mixing of orbitals with other j z values. 

In short, the present microscopic calculation demonstrates that apart from the two proton 
holes in the g$/2 orbital and one neutron particle in the hn/2 orbital, the other neutrons 
distributed over several subshells above the A = 50 core also contribute jointly to the AMR 
band in 105 Cd. The total angular momentum results from the alignment of the proton holes 
and the mixing within the neutron orbitals. Because of this strong mixing between the 

22 



neutron orbitals, the phenomenological interpretation in Ref. j^j without proper treatment 
of the core is not fully justified. 

E. Nucleon currents and densities 
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FIG. 10: (Color online) Dirac currents for 105 Cd in the yz plane perpendicular to the rotational 
axis at x = 1.473 fm for all neutrons (top), for the N = 50 core (middle), and for the seven valence 
neutrons (bottom) at Ml = 0.3 MeV (left panels) and Ml = 0.7 MeV (right panels). The shape 
and the size of the nucleus 105 Cd are represented by the solid lines corresponding to lines with 
constant neutron density p n = 0.01 fm -3 . 



As discussed in the introduction, the time reversal symmetry is broken in rotating nuclei. 
This leads to polarization currents and nuclear magnetism, i.e., time-odd components of the 
fields. In the following, we will investigate the current distributions. 
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FIG. 11: (Color online) Same as Fig. [10] but in the xz plane parallel to the rotational axis with 
y = 3.461 fm. 

In Fig. [10] we show Dirac currents for the neutrons in the nucleus 105 Cd. We present the 
currents in the yz plane at x = 1.473 fm for two different frequencies, in the left panels for 
HQ = 0.3 MeV and in the right panels for HQ = 0.7 MeV. The upper two panels (a) and (b) 
show the total currents of all the neutrons, the panels (c) and (d) in the middle show the 
neutron currents in the N = 50 core and the lower two panels (e) and (f) show the currents 
of the seven valence nucleons. 

It is clearly seen that the currents in the TV = 50 core are very small at Ml = 0.3 MeV and 
remain that even though the rotational frequency surge from HQ = 0.3 MeV to 0.7 MeV. 
Of course, the particles in the core feel also a strong Coriolis field OJ, which leads to 
considerable mixing of the different single particle wave functions. However, the core corre- 
sponds approximately to a closed spherical shell. It is represented by a Slater determinant 
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FIG. 12: (Color online) Dirac currents for 105 Cd in the yz plane with x = 1.473 fm for all protons 
(top), for the Z = 50 core (middle), and for the two valence proton holes (bottom) at Ml = 0.3 MeV 
(left panels) and Ml = 0.7 MeV (right panels). The shape and the size of the nucleus 105 Cd 
are represented by the solid lines which corresponds to to lines with constant proton density 
Pp = 0.01 fm" 3 . 

not depending on the mixing or the individual wave functions of the orbits. All the time- 
reversed partners are nearly equally occupied and this leads to a cancellation of the angular 
momenta and the currents. The net angular momenta and currents stay always negligible. 
In contrast, in the valence shell there are many empty orbits and the Coriolis field forces 
the seven valence neutrons to occupy only orbitals with angular momenta aligned in the 
^-direction. Orbits with opposite angular momenta are kept empty. Therefore, as shown 
in Fig. |H1 net alignments grow considerably with increasing angular velocity and the same 
holds for the currents. The nuclear magnetic fields induced by these currents enhance this 
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effect in a self-consistent way, because the isoscalar as well as the isovector current-current 
interactions between nucleons are always attractive. Repulsive current-current interactions 
between particles and anti-particles play no role because of the no-sea approximation. As 
a consequence, the total neutron currents are primarily generated by the valence neutrons 
and, as the nucleus is rotating around the x axis, the total current pattern of Fig. [TU] is 
dominated by a large central symmetric vortex in the yz plane. This current distribution is 
by no means irrotational. 
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FIG. 13: (Color online) Same as Fig. [12] but in the xz plane with y = 3.461 fm. 



In full analogy we show in Fig. [TTJ Dirac currents for 105 Cd in the xz plane parallel to 
the rotational axis at y — 3.461 fm. Contributions of all neutrons are given at the top, those 
of the N = 50 core in the middle, and those of the seven valence neutrons at the bottom. 
Two angular velocities Ml = 0.3 MeV (left panels) and HQ = 0.7 MeV (right panels) are 
presented. Again the total currents are dominated by the valence particles. However, since 
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the nucleus is rotating around the x axis, the current pattern in the xz plane is dominated 
by a longitude flow with negligible components in the x direction. 

In Fig. H21 same as in Fig. [TU] but for protons, the Dirac currents for 105 Cd are given 
for all protons (top), for the Z = 50 core (middle), and for the two valence proton holes 
(bottom) at htt = 0.3 MeV (left panels) and Ml = 0.7 MeV (right panels). 

Since there are only 48 protons in 105 Cd, the total currents are the difference between the 
contributions from Z = 50 core and the two valence proton holes. Similar to the currents 
of the N = 50 core, the currents for the Z = 50 core are also negligible at Ml =0.3 MeV or 
0.7 MeV due to the cancellation of the time reversal partners which also provides a negligible 
angular momentum (see Fig. [S]). As a consequence, the total proton currents are essentially 
build from the two valence holes which are enhanced with the growing frequency. However, 
since the total spin of protons is smaller than that of neutrons, the corresponding Coriolis 
term and the strengths of the proton currents are significantly weaker than the neutron ones 
at a given rotational frequency. 

In Fig. [131 same as in Fig. [12] but in the xz plane with y = 3.461 fm, the Dirac currents 
for 105 Cd are given for all protons (top), for the Z = 50 core (middle), and for the two 
valence proton holes (bottom) at Ml = 0.3 MeV (left panels) and Ml = 0.7 MeV (right 
panels). 

For the strengths, the currents in the xz plane increase with rotational frequency, which 
are similar as the currents in the yz plane. For the current pattern, similar as the neutron 
cases, the proton current pattern in the yz plane is dominated by a central symmetric vortex 
and in the xz plane by a appreciable longitude flow, which is connected with the fact that 
nucleus here is rotating around the x axis. However, different from the neutron currents in 
the xz plane, the x components of proton currents are no longer negligible. The difference 
can be understood from the so-called two-shears-like mechanism in Fig. [71 The two proton 
holes in orbital, which originally orient back to back along the z direction, gradually 
align towards the x axis and create flow not only in the longitude z direction but also in the 
x direction. 

It is interesting to investigate the evolution of the density distributions with increasing 
angular velocity and the relation between the currents and the density distributions. In 
Fig. [T41 we show the density distribution of 105 Cd in the yz plane at x — 1.473 fm for all 
neutrons at Ml = 0.3 MeV (panel a) and at Ml = 0.7 MeV (panel b). It is found that the 
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FIG. 14: (Color online) Density distribution for 105 Cd in the yz plane with x = 1.473 fm for all 
neutrons (top), N = 50 core (middle), and the seven valence neutrons (bottom) at Ml = 0.3 MeV 
(left panels) and Ml = 0.7 MeV (right panels). Note that different color scale is used for the valence 
neutrons for better visualization of densities. 

neutron density has a prolate-like shape which is symmetric around the z axis and changes 
with increasing rotational frequency quantitatively toward to a near-spherical distribution. 
This is consistent with the deformation evolutions shown in Fig. [6j 

The neutron density is composed of the N = 50 core as shown in panels (c) and (d) as 
well as the seven valence neutrons in the panels (e) and (f) in Fig. HH It is found that 
both the densities for the seven valence neutrons and the N = 50 core have prolate-like 
distribution with the symmetry axis z axis. Although the densities for the valence neutrons 
is one order of magnitude smaller than those of the N = 50 core, the neutron currents in 
Figs. [10] and [TTJ are mainly produced by the seven valence neutrons rather than the N = 50 
core. Examining the density distributions and the currents of the valence neutrons, it is 
interesting to note that the minima in density distributions correspond to the maxima of 
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FIG. 15: (Color online) Same as Fig. Q3]but in the xz plane with y = 3.461 fm. 
the current strengths. 

The same conclusion can also be reached by focusing on the neutron density distributions 
in the xz plane shown in Fig. [15], where the density distributions for all neutrons, N = 50 
core, and the seven valence neutrons in 105 Cd are displayed in the xz plane with y = 3.461 fm. 

In Fig. [16] same as in Fig. [14] but for protons, the density distributions for 105 Cd are 
displayed for all protons (top), for the Z = 50 core (middle), and the for the two valence 
proton holes (bottom) at Ml = 0.3 MeV (left panels) and Ml = 0.7 MeV (right panels). 

Similar as the neutron case, the proton density exhibits a prolate-like shape with z axis 
as the symmetry axis and evolutes toward a near-spherical distribution with the rotational 
frequency. Here the total proton density distribution can be decomposed into the density 
distributions of the Z = 50 core in panels (c) and (d) and the two valence proton holes in 
panels (e) and (f) in Fig. [16] It could be seen that the density distributions for the two 
valence proton holes have an oblate-like shape in contrast with the prolate shape for the 
Z = 50 core. By subtracting the density of two valence proton holes from the Z = 50 core, 
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FIG. 16: (Color online) Density distribution for 105 Cd in the yz plane with x = 1.473 fm for all 
protons (top), Z = 50 core (middle), and the two valence proton holes (bottom) at fr£i = 0.3 MeV 
(left panels) and Ml = 0.7 MeV (right panels). Note that different color scale is used for the valence 
proton holes for better visualization of densities. 

the resulting proton density is driven toward a more prolate shape. 

Again the densities for the valence proton holes are one order of magnitude smaller than 
those of the Z = 50 core, the proton currents in Figs. [12] and [13] are mainly produced by the 
valence proton holes rather than the Z = 50 core. Examining the density distributions and 
the currents of the valence proton holes, it is interesting to note that the maxima in density 
distributions correspond to the maxima of the current strengths due to the hole character. 

The same arguments also apply for the proton density distributions shown in Fig. [TTl 
where the density distributions of all protons, Z = 50 core, and the two valence proton holes 
in 105 Cd are displayed in the xz plane with y = 3.461 fm. 
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FIG. 17: (Color online) Same as Fig. [16] but in the xz plane with y = 3.461 fm. 
V. SUMMARY 

In summary, covariant density functional theory is used for a systematic investigation 
and a detailed analysis of the AMR band in the nucleus 105 Cd. Following a short report [l| 
these calculations within the framework of the TAC model present the first fully microscopic 
and self-consistent description of a AMR band in nuclei. 

Using the point-coupling density functional PC-PK1, the configuration for the observed 
AMR band in 105 Cd is fixed by analyzing the single-particle Routhians. With the con- 
figuration thus obtained, the tilt angle 6q for a given rotational frequency is determined 
self-consistently by minimizing the total Routhian with respect to the angle 6q. It demon- 
strates in an unambiguous way that the rotational axis is always along the x axis. 

The energy spectrum, the total angular momentum, the kinetic and dynamic moments 
of inertia, and the B(E2) values of the AMR band in 105 Cd are calculated and agreement is 
found with available experimental data. The change in the energy along the spectrum in the 
band comes to a large part from the kinetic energy but with a moderate competition from 
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the interaction energy. The observed increase of the 3^/B(E2) ratios is well reproduced 
reflecting the fact that the moment of inertia 3^ is essentially constant whereas the B(E2) 
values rapidly decline along the band. This fact clearly distinguishes AMR bands from 
smoothly terminating bands. 

In contrast to the phenomenological model in Ref. [u]] where the nucleus is divided artifi- 
cially into a core and serval valence particles or holes, in the present microscopic calculation 
both the energy and the angular momentum along a band come entirely from individual 
nucleons. By investigating microscopically the contributions of neutron and protons to the 
total angular momentum, the two-shears-like mechanism in the AMR band is clearly illus- 
trated. 

Finally, the Dirac currents are presented and discussed in details for the AMR band in 
105 Cd. They are originally induced by the Coriolis operator, but considerably increased 
by the nuclear magnetic fields resulting from a self-consistent polarization mechanism. By 
decomposing the currents into the contributions from the core and from the valence particles 
or holes, it is found that the core stays mostly inert. The structure of the total currents is 
essentially determined by the valence particles or holes. Their size and their spatial pattern 
depends on the specific single-particle orbitals and, of course, on the rotational frequency. 
In addition we discuss the correlations between the currents and the density distributions 
in the present calculations. 
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